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ABSTRACT 

This  report  outlines  accuracy  benefits  in  the  time  difference  of  arrival  (TDOA) 
estimation  of  emitters  employing  2-FSK  signals  over  a  corresponding  OOK  sig¬ 
nal.  We  consider  TDOA  estimation  via  the  conventional  approach,  using  max¬ 
imization  of  the  magnitude  of  the  complex  correlation  function.  By  drawing  an 
analogy  with  the  improvement  that  is  theoretically  possible  by  incorporating 
phase  information  into  TDOA  estimation  for  an  arbitrary  signal,  it  is  suggested 
that  for  sufficiently  high  SNR,  substantially  more  accurate  estimation  of  the 
TDOA  is  possible  with  an  FSK  signal. 
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Considerations  on  the  Relative  Accuracy  of  TDOA 
Estimation  for  FSK  and  OOK  Signalling 

Executive  Summary 

The  estimation  of  the  relative  time  delay  between  a  signal  component  common  to  that 
received  at  two  spatially  separated  antennas  is  a  problem  of  common  interest.  The  time 
difference  of  arrival  (TDOA)  is  usually  estimated  as  the  argument  of  the  peak  in  the  cross¬ 
correlation  computed  between  the  two  received  signals.  Conventionally,  to  minimize  the 
sampling  rate  and  consequently  the  computational  complexity,  the  correlation  is  performed 
on  the  complex  envelopes  of  the  received  signals. 

This  report  compares  the  TDOA  estimation  variance  (accuracy)  achievable  with  two 
popular  modulation  schemes:  on-off  keying  (OOK),  and  2-mode  frequency  shift  keying  (2- 
FSK).  The  former  is  characterised  by  the  carrier  frequency,  fc,  and  the  bit  rate  Rb]  while 
the  latter  is  also  characterised  by  the  frequency  deviation  fo  (of  the  two  FSK  components 
from  fc). 

The  Cramer- Rao  bound  for  the  variance  of  the  TDOA  is  inversely  dependent  on  an 
integral  (over  frequency)  involving  the  ‘SNR  spectra’,  the  cross-correlation  integration 
time,  the  signal  bandwidth  as  well  as  a  quadratic  frequency  multiplier.  For  the  case  in 
which  the  signals  and  noise  have  (piecewise)  flat  power  spectra,  the  latter  term  leads  to 
the  TDOA  estimation  variance  being  inversely  proportional  to  the  Tms  radian  frequency’. 
This  suggests  that  the  2-FSK  signal  type  should  have  a  relatively  lower  TDOA  estimation 
variance  than  the  OOK  signal  type,  with  the  relative  improvement  growing  with  fo/Rb- 
This  is  reflected  in  the  complex  correlation  envelope  by  the  2-FSK  signal  showing  a  rela¬ 
tively  sharper  primary  correlation  peak,  by  a  factor  related  to  fo/Rb-  On  the  other  hand, 
the  complex  correlation  envelope  of  the  2-FSK  signal  includes  an  oscillatory  component, 
leading  to  ambiguous  secondary  peaks,  lying  at  multiples  of  l/2/o  from  the  true  TDOA. 
Heuristically,  the  likelihood  of  such  secondary  peaks  causing  erroneous  TDOA  estimation 
increases  within  increasing  fo/Rb  and  decreasing  SNR. 

With  the  aim  of  obtaining  a  stronger  understanding  relating  to  SNR,  fo/Rb  and  the 
relative  improvement  in  TDOA  estimation  provided  by  2-FSK  over  OOK,  we  draw  on  an 
analogy  with  the  TDOA  estimation  improvements  achievable  by  incorporating  complex 
correlation  phase  information  with  the  complex  correlation  envelope.  A  series  of  theoretical 
studies  have  been  published  on  this  analogous  consideration.  These  imply  that,  for  the 
2-FSK  vs  OOK  comparison  case,  two  SNR  thresholds  SNR2  >  SNR\  exist.  In  particular: 
(i)  for  SNR  <  SNRi,  the  likelihood  of  ambiguity  induced  TDOA  estimation  errors  for 
the  2-FSK  signal  type  is  substantial,  leading  to  essentially  the  same  TDOA  estimation 
accuracies  as  that  provided  by  OOK;  (ii)  for  SNR  >  SNR-2,  the  likelihood  of  ambiguity 
induced  TDOA  estimation  errors  with  the  2-FSK  signal  type  is  negligible,  leading  to 
this  signal  type  showing  substantially  improved  TDOA  estimation  over  the  OOK  signal 
type;  (iii)  the  SNR  thresholds  SNR\  and  SNR2  increase  as  fo/Rb  increases.  The  report 
concludes  with  a  simulation  study  confirming  the  above  features.  Importantly,  in  many 
applications,  the  required  SNR  (>  SNR2)  to  avoid  ambiguity  induced  errors  is  practically 
achievable,  leading  to  the  2-FSK  signal  type  showing  an  improvement  (reduction)  in  TDOA 
estimation  variance  on  the  order  of  (fo/Rb)2- 
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1  Introduction 


The  estimation  of  the  relative  delay  between  a  signal  component  common  to  that  received 
at  two  spatially  separated  antennas  is  a  problem  of  general  interest.  Consider  the  reception 
of  a  transmitted  signal  u(t)  by  two  spatially  separated  antennas,  denoted  r\(t)  and  r2{t) 

ri(t)  =  a\ u(t  -  Ti )  +  n±(t) 

r2{t)  =  a2u(t  -  T2)  +  n2(t),  (1) 

where  ni(t)  and  n2(t)  are  additive  noise  with  arbitrary  distribution.  For  clarity  and  to 
simplify  the  analysis,  the  scalar  coefficients  a\  and  a2  will  arbitrarily  be  set  to  one  and 
T\  will  be  set  equal  to  zero,  so  that  the  relative  time  delay,  sometimes  known  as  the  Time 
Difference  of  Arrival  (TDOA)  or  Differential  Time  Offset  (DTO)  [3],  T  =  T2.  This  is  often 
estimated  as  the  argument  of  the  correlation  peak  T  =  argmaxr  Rrir2(r).  If  one  had 
the  ability  to  exactly  compute  this  cross  correlation  function  then  the  approach  would  be 
optimal,  but  when  the  correlation  must  be  estimated,  the  optimal  (maximum  likelihood) 
approach  is  to  emphasise  the  energy  associated  with  the  signal  component  by  appropriately 
pre-filtering  the  r\(t)  and  r2{t)  in  a  generalised  correlation  [1,  5].  In  the  absence  of  a  priori 
knowledge  of  the  signal  and  noise  characteristics,  a  pragmatic  approach  can  be  as  simple 
as  prefiltering  the  signals  to  approximately  the  bandwidth  that  contains  most  of  the  signal 
energy. 

To  minimise  the  sampling  rate  and  consequently  the  computational  complexity,  cor¬ 
relation  is  often  performed  [5]  on  the  complex  envelopes  fi(t)  and  r2{t)  so  that  T  = 
argmaxr  |Ar1f2('r)|-  Consider  the  case  where  the  signal  and  noise  are  restricted  to  within 
a  passband  bandwidth  fc±B/ 2,  where  fc  and  B  are  the  carrier  frequency  and  bandwidth 
respectively  and  where  the  correlation  is  estimated  over  an  L  second  interval.  For  suffi¬ 
ciently  high  signal  to  noise  ratio  (SNR)1  the  Cramer  Rao  bound  for  the  variance  of  the 
estimation  of  time  delay  using  complex  envelopes  is  given  by  [5] 


a 


2 

min 


rfc+B/2 

4tt2L  /  2 SNR(f)(f  -  fc)2df 

■Jfc-B/2 


-1 


where 


SNR(f) 


Guu(f)/Gnlnl(f).Guu(f)/Gn2n2(f) 

1  +  Guu(f)/Gnlnl(f)  +  Guu(f)/Gn2n2(f) 


(2) 


(3) 


and  Guu(f),  Gnini{f )  and  Gn2n2(f)  are  the  signal  and  noise  auto-spectra.  It  can  be  shown 
[2]  that  for  the  case  where  both  the  signal  and  noise  have  white  power  spectra,  the  variance 
is  bounded  by  the  minimum  [3] 


(T 


2 

min 


l 

Wills' 


(4) 


lrTliere  exists  an  SNR  threshold  region  [4],  below  which  time  delay  estimation  contains  ambiguities,  but 
above  which,  the  variance  is  bounded  by  the  Cramer  Rao  low  bound. 
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where  the  equivalent  SNR 


1  1  /  1  1  1  \ 

—  —  -  ( - 1 - 1 - )  (5) 

leq  2  72  7172/ 

embodies  the  SNR  associated  with  the  signal  received  at  the  two  apertures,  71  and  72.  (3 
is  defined  as  the  ‘mis  radian  frequency’  [3], 

V2 

,  (6) 


(3  =  2n 

where  the  Guu(f )  =  Guu(f  -  /c). 


fZcfGuuUW 
IZ  Guu(f)df 


Although  equation  (4)  applies  to  the  white  signal  and  noise  case,  it  nevertheless  sug¬ 
gests  more  general  dependencies  for  the  achievable  TDOA  estimation  variance.  In  partic¬ 
ular,  it  can  be  seen  that  the  variance  will  reduce  as  the  bandwidth  of  the  signal  increases, 
as  might  be  expected.  However,  it  doesn’t  specifically  reveal  whether  there  are  benefits 
associated  with  any  particular  modulation  format.  In  the  next  section,  an  analysis  will  be 
presented  that  makes  clear  the  variance  reduction  that  could  be  achieved  by  conducting 
the  correlation  estimate  at  an  RF  frequency.  Section  3  leverages  from  that  analysis  to 
suggest  that  FSK  modulation  exhibits  a  property  that  provides  a  TDOA  estimate  that 
can  be  more  accurate  than  a  comparable  OOK  modulation.  The  accuracy  improvement 
is  illustrated  through  simulations. 


2  Passband  TDOA  estimation 

Although  it  is  common  for  TDOA  to  be  estimated  from  the  computation  of  the  cross 
correlation  between  two  complex  baseband  signals,  it  is  theoretically  possible  to  obtain  an 
improved  estimate  by  instead  computing  the  cross  correlation  between  the  signals  at  their 
received  RF  frequencies.  At  sufficiently  high  SNR,  the  estimation  variance  is  governed  by 
the  Cramer  Rao  bound,  which  for  a  real,  passband  signal  with  fc  Hz  carrier  frequency,  is 
given  by  [5] 

-1 

(7) 

In  comparing  (7)  to  (2),  it  can  be  seen  that  the  bound  differs  only  by  the  inclusion  of  f2 
in  lieu  of  (/  —  /c)2.  For  complex  processing,  the  variance  depends  only  on  the  frequency 
spread  of  the  signal  about  its  carrier  frequency,  whereas  for  real  processing,  the  variance 
is  reduced  as  the  SOI  frequency  content  increases.  In  sections  2.1  and  2.2  an  analysis  is 
presented  that  provides  insight  into  the  reason  for  this. 


rfc+B/2 

47 T2L  /  2 SNR(f)f2df 

Jfc-B/2 


2.1  Complex  baseband  TDOA  estimation 


Consider  the  bandpass  representation  of  the  signal 

u(t)  =  3?  [s(t)  exp(j27r/ct)] 

=  a(t)  cos(27r/ct)  —  b(t)  sin(27r fct), 


(8) 
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where  s(t)  =  a(t)+jb(t)  and  which  is  assumed  to  be  transmitted  from  a  fixed  location.  For 
the  ensuing  analysis,  we  neglect  the  noise  components  on  the  received  signals,  mindful  that 
these  will  give  rise  to  terms  that  only  distract  from  the  observations.  Complex  baseband 
TDOA  estimation  considers  the  complex  quasi-basebanded,  downconverted  signals: 

n(t) — *  xi(t)  =  s(t)exp(j2nfIt)  (9) 

r2 (t) — >  x2(t)  =  s(t  —  T)  exp(j'27r//f)  exp(— j27r/cT)  (10) 


where  fj  is  the  quasi-basebanded  intermediate  frequency.  Estimation  of  the  TDOA  is  via 
maximisation  of  the  magnitude  of  the  correlation  function: 


T  = 


arg  max  |  R(r)  \ 

T 

=  argmaxl  /  x\{t  —  r)x2(t)dt\ 

T  Jt 

=  arg  max  |  exp(j27r(/cT  -  far))  J s(t  -  r)s*(t  -  T)dt\ 
=  arg  max  |  J  s(t  —  r)s*(t  —  T)dt\ 


(11) 


As  an  example,  consider  an  emitter  employing  OOK  modulation  such  that  s(t)  =  a(t)  is 
a  binary  unipolar  NRZ  signal,  with  amplitude  A  and  bit  rate  Rb-  For  simplicity,  we  set 
A  =  1.  Then  the  TDOA  estimate  is  from  (11), 


f 


arg  max 

T 


r)a(t  —  T)dt 


(12) 


2.2  Real  Passband  TDOA  Estimation 


The  complex  baseband  approach  from  section  2.1  ignores  the  phase  information  of  R{r). 
One  approach  to  exploiting  the  phase  information  is  to  process  the  real  passband  collected 
signals  which  again,  neglecting  the  noise  terms  can  be  represented  as 

r\(t)  =  a(t)  cos(2ir fct)  —  b(t)  sin(27r/cf)  (13) 

r2{t)  =  a(t  -  T )  cos(27r/c(t  -  T))  -  b(t  -  T )  sin(27r/c(t  -  T)).  (14) 

Estimation  of  the  TDOA  is  then  via 


T  = 


arg  max  R(t) 

T 

=  arg  max  /  r\(t  —  T)r2(t)dt 

T  Jt 

=  arg  max[cos(27r fc(T  —  r))  ~  x)a(t  —  T)  +  b{t  —  r)b{t  —  T)dt 

+  sin(27r fc(T  —  r))  ~  T)b{J  —  T)  —  b(t  —  r)a(t  —  T)dt], 


where  terms  at  double  the  carrier  frequency  have  been  assumed  to  integrate  out.  It  is 
again  useful  to  consider  the  OOK  example  from  section  2.1.  In  this  case,  the  TDOA 
estimate  simplifies  to 


f 


arg  max 

T 


cos(27 t/c(T 


r)) 


r)a{t  —  T ) 


dt. 


(15) 
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The  additional  “phase  related”  terms  cos(27 r/c(T  — r))  and  sin(27r /c(T  —  r))  potentially 
lead  to  significantly  improved  TDOA  estimation,  due  to  a  modulation  sharpening  effect 
on  the  peak  of  the  correlation  function.  The  improvement  is  dependent  on  the  ratio  of 
fc/B.  These  “phase”  terms  however,  also  introduce  estimation  ambiguities: 

T  =  T-  n/fc  (16) 

where  n  is  an  integer.  The  susceptibility  to  these  ambiguities  increases  as  fc/B  increases 
and  as  the  SNR  and  the  correlation  integration  time  (L)  decreases. 

As  an  illustration,  consider  the  case  where  the  signal  and  noise  are  both  white  within 
fc  ±  B/2.  For  sufficiently  strong  SNR,  satisfying  [4,  5] 


SNR 

>  6f^/(TT2B3L)[(j)~1(B2/ (24/c))]2 

(17) 

where  4>(y) 

/*oo 

=  1/V2t:  J  exp(—x2 /2)dx, 

(18) 

the  improvement  (reduction)  in  TDOA  estimation  variance  (using  the  real  passband  pro¬ 
cessing  approach  rather  than  the  conventional  complex  baseband  processing  approach),  is 
given  by  [5] 

aT,conv/aT  =  1  +  12  {fc/B)2.  (19) 

Although  not  exactly  applicable  to  more  general  spectral  shapes,  these  expressions  never¬ 
theless  suggest  the  order  of  improvement  that  can  be  obtained  by  estimating  the  TDOA 
using  the  passband  approach.  It  should  be  noted  however,  that  this  improvement  assumes 
operation  in  the  SNR  region  governed  by  the  Cramer  Rao  bound  and  that  any  ambiguity 
induced  TDOA  estimates  are  clearly  detected  as  outliers  and  are  ignored.  As  the  SNR 
falls  (e.g.  for  the  case  of  bandlimited  white  noise  and  signal,  below  the  threshold  of  (17)), 
the  relative  improvement  decreases  due  to  the  increased  influence  of  ambiguity  induced 
TDOA  errors. 

For  sufficiently  poor  SNR  (e.g.  for  the  case  of  bandlimited  white  noise  and  signal), 

SNR  <  2.76fc/(TT2B3L)  (20) 

the  influence  of  the  TDOA  ambiguities  dominates  and  the  two  approaches  lead  to  essen¬ 
tially  the  same  TDOA  estimation  variance. 

To  re-iterate,  the  real  passband  processing  approach  provides,  at  worst  (corresponding 
to  ‘poor’  SNR  conditions  as  described  by  equation  (20),  TDOA  estimation  accuracies  equal 
to  that  provided  by  the  conventional  complex  baseband  processing  approach.  For  suffi¬ 
ciently  ‘good’  SNR  conditions,  the  real  passband  processing  approach  provides  improved 
accuracy. 

Simulations  were  performed  to  conduct  a  spot  check  on  the  validity  of  equation  (19) 
and  the  extent  of  its  applicability  beyond  a  bandlimited  white  signal  to  a  modulated 
waveform  such  as  OOK.  For  this  purpose,  the  following  parameters  were  considered:  an 
SNR  equal  to  -10  dB  within  a  B  =  10kHz  signal  bandwidth,  a  1.5  sec  integration  period 
and  fs  =  312.5  kHz.  With  a  carrier  frequency  fc  =  5  kHz,  the  respective  mis  TDOA  for 
real  passband  and  complex  TDOA  estimation  was  1.9//sec  and  3.2  fisec.  The  improvement 
in  simulation  variance  is  then  approximately  a  ratio  of  2.8,  whereas  equation  (19)  would 
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suggest  a  ratio  equal  to  4.0.  At  fc  =  10  kHz  the  mis  TDOA  for  real  passband  estimation 
was  0.9//sec  while  rms  TDOA  for  complex  estimation  was  again  3.2  /rsec.  In  this  case,  the 
variance  ratio  was  approximately  12.6,  closely  agreeing  with  13.0  as  suggested  by  (19). 
These  results,  although  by  no  means  conclusive,  give  confidence  to  the  validity  of  equation 
(19). 


Now  consider  simulations  using  an  OOK  modulated  signal,  bandlimited  to  its  baud 
rate  frequency  fb  =  5  kHz,  so  that  B  =  10  kHz  again.  Under  the  same  conditions  as 
the  previous  simulations  for  the  white  signal,  with  fc= 5  kHz,  the  rms  TDOA  for  the  real 
passband  and  complex  estimation  was  2.4  [i sec  and  14.3  ^xsec  respectively2.  The  variance 
ratio  was  consequently  approximately  35.5,  substantially  greater  than  for  the  bandlimited 
white  signal.  With  fc  =  10kHz,  the  rms  TDOA  for  real  and  complex  estimation  was  0.96 
and  14.3  jusec  respectively,  with  a  variance  ratio  approximately  equal  to  222.  Clearly,  the 
degree  of  improvement  associated  with  an  OOK  modulated  signal  is  substantially  greater 
than  with  a  bandlimited  white  signal. 


The  reason  that  the  degree  of  TDOA  accuracy  improvement  obtained  through  esti¬ 
mating  the  correlation  function  using  the  real,  passband  data  is  substantially  better  for 
a  bandlimited  OOK  signal  than  for  a  comparable  bandlimited  white  noise  signal  relate 
specifically  to  the  shape  of  the  correlation  function  of  the  message  in  each  case.  In  par¬ 
ticular,  the  correlation  function  of  a  bandlimited  square  waveform  (the  message  of  an 
OOK  signal)  is  broader  than  for  a  bandlimited  white  noise  waveform  3.  Because  of  this 
greater  broadness,  the  sharpening  effected  by  the  modulation  has  more  impact  than  for 
the  bandlimited  white  noise  signal. 


In  Appendix  A  it  is  shown  that  the  improvement  obtained  through  computing  the 
correlation  function  from  the  passband  data  can  be  alternatively  obtained  through  a  more 
explicit  incorporation  of  phase  information  in  a  modified  complex  baseband  approach. 
It  needs  to  be  emphasised  that  either  mechanism  to  achieve  variance  reduction  through 
the  exploitation  of  phase  information  requires  a  sampling  rate  fs  >  2  fc  which,  typically, 
is  significantly  higher  than  the  sampling  rate,  fs  >  2 B,  required  with  the  conventional 
complex  baseband  approach. 


2Note  that  at  this  SNR  and  with  these  parameters,  performance  is  at  the  edge  of  the  region  bounded 
by  the  Cramer  Rao  bound  and  anomalous  correlation  peaks  start  to  be  observed  in  the  simulations. 
Consequently,  TDOA  estimates  were  hard  limited  to  between  ±10  samples  from  the  true  peak,  to  ensure 
the  results  can  be  interpreted  as  applicable  to  the  Cramer  Rao  bounded  region. 

3  Specifically,  first  consider  the  OOK  signal.  The  square  message  has  a  power  spectrum  defined  by  a 
sine  squared  shape,  with  first  nulls  at  ±13/2.  Its  autocorrelation  function  is  the  inverse  Fourier  transform, 
which  is  a  triangular  function  that  tapers  to  zero  at  r  =  ±2/13.  The  bandlimiting  of  the  OOK  signal 
constitutes  multiplication  of  its  power  spectrum  with  a  square  filter,  extending  to  ±13/2  Hz,  which  has  an 
inverse  Fourier  transform  that  is  a  sine  function  with  first  nulls  at  ±1/13.  Consequently,  the  autocorrelation 
function  of  the  bandlimited  OOK  message  is  spread  by  the  convolution  of  the  sine  and  triangular  function. 
The  bandlimited  white  noise  signal,  on  the  other  hand,  has  an  autocorrelation  function  that  is  simply  that 
same  sine  function. 
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3  Comparison  between  TDOA  estimation 
accuracy  for  OOK  and  FSK  modulation 


In  section  2.2  it  was  shown  that  for  an  arbitrary  signal,  if  the  TDOA  is  estimated  using 
a  passband  correlation,  the  resulting  correlation  function  can  in  some  circumstances  be 
sharpened  by  an  inherent  oscillatory  term.  In  this  section,  it  will  be  shown  that  for  the 
specific  case  of  an  FSK  signal,  an  analogous  sharpening  can  also  be  achieved,  relative  to 
an  OOK  signal,  even  when  using  a  conventional  complex  baseband  correlation  estimate. 

Consider  an  emitter  employing  a  2-FSK  signal  formed  as  a  composite  of  two  phase- 
inverted  OOK  signals: 

u(t)  =  m(t)  cos(27t(/c  +  f0)t)  +  (1  —  m(t))  cos(27r(/c  —  fQ)t  —  4>) 

=  m(t)  cos(27t fat)  cos(27r/cf)  —  m(t)  sin(27r/0f)  sin(27r/cf) 

+(1  —  m(t))  cos(27r/Df  +  (j>)  cos(2irfct)  +  (1  —  m(t))  sin(27r/0t  +  0)  sin(27r fct) 

=  Re[s(t)  exp(j27r/cf)] 

where  4>  is  an  arbitrary  phase  offset,  m(t)  is  a  binary,  unipolar,  NRZ  message  and 

s(t)  =  m(t )  exp(j27r/0f)  +  (1  -  m(t))  exp(-j27r/0f  -  j<f>).  (21) 


The  signal  is  again  collected  by  two  stationary  apertures,  as  in  equation  (1).  Complex 
baseband  estimation  of  the  TDOA,  T  is  achieved  via  the  downconverted  signals  which, 
again  ignoring  noise  contributions  are 

ri(t)  — »  xi (t)  =  s(t)  exp(j'27r/jt) 

r2(t)  — >  x2 (t)  =  s(t  -  T)  exp(j27r fit)  exp(— j27r/cT).  (22) 


The  TDOA  is  estimated  using 

T  =  arg  max  |  R(t)  \ 

T 


=  arg  max 

r 

=  arg  max 


xi(t  —  r)x2(t)(it| 
exp(j2vr(/cT- //r))  f  s(t 


r)s*(t  —  T)dt 


(23) 


The  exponential  term  in  (23)  has  no  effect  on  the  magnitude,  so  can  be  omitted.  Further, 
the  integrand  can  be  expanded  as 


s(t  -  r)s*(f  -  T)  =  ( m(t  -  r)  exp(j2nf0(t  -  r))  +  (1  -  m(t  -  r))  exp(-j2vr/0(t  -  r)  -  j<j>)) 

x  (m(t  -  T)  exp(j27r/0(t  -  T))  +  (1  -  m(t  -  T))  exp(-j2nf0(t  -  T)  -  j(f> ))* 
=  (m(t  -  t)  exp(j27r/0(f  -  r))  +  (1  -  m(t  -  r))  exp(-j2vr/0(t  -  r)  -  j(j>)) 
x  (m(t  -  T)  exp(-j2vr/0(t  -  T))  +  (1  -  m(t  -  T))  exp(j27r/0(t  -  T)  +  j(j>)) 
=  m(t  —  T)m(t  —  T)  exp(j2irf0(T  —  r)) 

+(1  -  m(t  -  r))(l  -  m(t  -  T))  exp(-j27r/0(T  -  r)) 

+m{t  —  t)(1  —  m(t  —  T))  exp(j27r2f0t)  exp(— j2ir(T  +  r)  +  jty) 

+m(t  -T)(  1  -  m(t  -  t))  exp(-j27r2 fat)  exp(j27r(T  +  r)  -  j<j>) 
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If  it  is  assumed  that  the  terms  corresponding  to  double  the  FSK  frequency  deviation 
integrate  out,  then  the  integral  simplifies  to 


J s(t  —  r)s*(t  —  T)dt  =  J m(t  —  T)m(t  —  T )  exp(j27r/0(T  —  r)) 

+(1  —  m{t  —  r))(l  —  m(t  —  T))  exp(— j'27t/0(T  —  r))dt 
=  J m(t  -  t )m(t  -  T)  (exp(j2vr/0(T  -  r))  +  exp(-j27r/0(T  -  r))) 
+(1  —  m.(t  —  t)  —  m{t  —  T ))  exp(— j27r/0(T  —  r))dt 
=  2  cos(27t/o(T  —  t))  J m(t  —  r)m(t  —  T ) 

+  exp(—  j2irf0(T  —  r))  J (1  —  m(t  —  r)  —  m(t  —  T))dt. 


If  m{t)  has  a  balanced  distribution  of  logical  ‘l’s  and  ‘0’s  then  the  integration  of  the  term 
1  —  m{t  —  t)  —  m{  1  —  T)  approaches  zero  for  L»T,r.  Consequently,  the  integral  simplifies 
further  to 


J  s(t  —  r)s*(t  —  T)dt 


and  the  TDOA  estimate  is 


2  cos(27t  f0(T 


t)) 


r)m(t  —  T)dt 


f 


arg  max 

T 


^2cos(27 r/0(T 


t)) 


T)m{t  —  T)dt 


(24) 


This  has  a  structure  that  is  very  similar  to  equation  (15)  for  generic  real  passband  TDOA 
estimation,  with  the  correlation  magnitude  |-R(t)|  gain  maximized  by  T  =  r  =  T.  In 
an  analogous  manner  to  the  generic  passband  processing,  the  presence  of  the  cos(.)  term 
increases  the  curvature  (or  “sharpness”)  at  the  peak  value  (at  T  =  r  =  T).  which  can 
lead  to  enhanced  TDOA  estimation  accuracy.  Instead  of  being  dependent  on  fc,  for  the 
present  case  the  curvature  will  depend  on  /o .  The  magnitude  of  this  increase  in  curvature 
will  then  depend  on  the  factor  f0/B,  where  B  is  a  measure  of  the  signal  bandwidth  which, 
drawing  from  the  analogy  between  (15)  and  (24)  is  suggested  as  B  =  2 fo.  4. 

The  discussions  of  the  preceding  paragraph  have  exploited  the  similarities  between 
(15)  and  (24)  and  have  ignored  the  differences.  One  such  difference  is  the  presence  of  the 
factor  ‘2’  in  the  expression,  which  represents  a  3dB  greater  correlation  peak  than  the  OOK 
peak  correlation  value,  but  it  should  be  appreciated  that  this  is  due  to  the  FSK  signal 
having  double  the  power.  The  other  significant  difference  in  the  expressions  is  the  presence 
of  the  magnitude  operation  in  (24).  One  of  the  clearest  consequences  of  the  magnitude 
operation  will  be  a  dependence  on  the  double  frequency  term,  2/o,  due  to  its  rectifying 
action.  Consequently,  the  oscillatory  nature  of  the  cos(.)  term  will  introduce  potential 
TDOA  estimation  ambiguities  at 


T  =  T±  n/2f0 

4The  analysis  for  the  generic  passband  TDOA  estimation  assumes  the  signal  and  noise  are  both  white 
within  the  band  fc  ±  B/ 2.  An  exact  analogy  between  that  analysis  and  the  present  would  require  white 
signal  and  noise  between  /o  ±  B / 2. 


UNCLASSIFIED 


7 


DSTO-TR-2516 


UNCLASSIFIED 


where  n  is  an  integer.  The  susceptibility  to  these  ambiguities  increases  as  f0/B  increases 
and  as  the  SNR  and  the  correlation  integration  time  (L)  decreases.  However,  unlike  for  the 
generic  passband  case,  the  sharpening  that  accompanies  the  FSK  signal  is  dependent  on  the 
frequency  deviation  and  not  the  RF  carrier  frequency.  Consequences  of  this  relationship 
include  it  being  more  likely  that  the  performance  improvement  is  realisable  in  practise,  as 
well  as  having  substantially  less  demanding  sampling  rate  requirements. 

By  drawing  on  the  similarities  between  equations  (15)  and  (24)5,  and  with  analogy  to 
(17),  a  rule  of  thumb  could  be  suggested  that  for  sufficiently  strong  SNR,  satisfying: 


SNR 

>  6fg/(TT2B3L)[(j)~1(B2/ (24/q  ))]2 

(25) 

where  4>(y) 

roo 

=  \/y/2n  J  exp(— x2 /2)dx 

(26) 

the  2-FSK  signal  would  provide  an  improved  (reduced)  TDOA  estimation  variance  (over 
the  OOK  signal)  given  by: 


aT,OOK/aT,2FSK  —  1  +  1 2(/0/F>)2.  (27) 

It  needs  to  be  acknowledged  that  the  OOK  and  FSK  signals  will  generally  have  power 
spectra  that  differ  significantly  from  being  white,  so  it  is  expected  that  it  will  not  be 
possible  to  apply  equation  (27)  exactly  in  practice.  It  should  also  be  noted  that  this 
improvement  assumes  that  any  ambiguity  induced  TDOA  estimates  are  clearly  detected 
as  outliers  and  are  ignored.  As  the  SNR  falls  below  the  threshold  of  (25),  the  relative 
improvement  decreases  due  to  the  increased  likelihood  that  not  all  ambiguity  induced 
TDOA  estimates  will  be  ignored. 

For  sufficiently  small  SNR,  the  two  signal  types  lead  to  essentially  the  same  TDOA 
estimation  variance.  With  analogy  to  equation  (20)  for  the  real  passband  TDOA  process, 
a  rule  of  thumb  could  be  suggested  that  the  2-FSK  signal  will  lead  to  improved  TDOA 
estimation  accuracy  if 

SNR  >  2.76/o/(7t2R3L).  (28) 

3.1  Illustration  through  simulation  results 

The  observations  of  the  previous  sections  are  illustrated  here  through  Matlab  simulation 
results.  Unless  otherwise  stated,  the  results  were  generated  with  with  arbitrarily  chosen 
parameters:  fs  =  312.5  kHz,  /&  =  5  kHz,  /o  =  10  kHz,  L  =  1.5  sec  and  a  carrier  frequency 
randomly  chosen  within  the  range  10  Hz  to  333  Hz.  The  bandwidths  of  the  signals  and 
noise  were  constrained  to  approximately  ±/&  Hz.  The  TDOA  was  randomly  selected  within 
a  ±1000  sample  range  and  100  iterations  were  used  in  each  trial.  The  SNR  is  measured 
within  the  signal  bandwidth.  That  is,  approximately  fc  ±  /),  for  OOK  and  a  2/&  band 
centered  at  both  fc  ±  /o  for  FSK. 

Substituting  fo  instead  of  /c.  It  could  be  conjectured  that  2/o  would  be  more  appropriate,  given  the 
rectification  effect  of  the  magnitude  operator  in  (24).  However  the  rectification,  although  giving  rise  to 
double  frequency  components,  will  not  change  the  shape  of  the  modulating  cosine  term  associated  with 
the  main  correlation  lobe. 
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Figures  1  and  2  show  a  central  truncation  of  the  correlation  functions  computed  for 
a  particular,  arbitrary  data  set,  for  OOK  and  FSK  signals  with  both  5  kHz  and  20  kHz 
frequency  deviation.  The  figures  correspond  respectively  to  a  received  SNR  of  +60dB  and 
-9dB,  on  both  signals.  The  sharpening  of  the  main  correlation  peak  for  FSK  with  respect 
to  OOK  is  clear.  However,  particularly  for  the  -9dB  case,  it  is  also  clear  that  there  is  an 
ambiguity  associated  with  the  selection  of  the  ‘correct’  correlation  peak. 

The  relative  mis  TDOA  estimation  accuracy  between  OOK  and  FSK  is  shown  in 
Figures  3  and  4.  Figure  3  shows  the  performance  over  the  approximate  received  SNR 
range  -14  dB  to  0  dB  and  shows  a  considerable  accuracy  advantage  with  FSK.  Previous 
experience  with  other  OOK  simulations  suggest  that  for  this  parameter  set,  the  detection 
of  anomalous  TDOA  peaks  begins  to  become  noticeable  below  -12  dB6,  so  this  represents 
the  lower  SNR  end  of  useful  TDOA  estimation. 

At  higher  SNR  levels,  smaller  deviations  are  associated  with  the  TDOA  estimation, 
so  a  faster  sampling  rate  is  required  to  resolve  such  performance.  Figure  4  shows  the 
simulation  results  with  fs  =  10  MHz.  In  order  to  achieve  computationally  tractable 
Matlab  simulations,  the  integration  period  was  shortened  to  L  =  0.15  sec  in  association 
with  the  higher  sampling  rate.  This  reduces  the  post-integration  effective  SNR,  so  that 
the  estimation  variance  at  0  dB  SNR  for  L  =  1.5  sec  should  roughly  equate  to  that  at 
8.5  dB  SNR  with  L  =  0.15  sec7.  Another  practicality  with  the  use  of  fs  =  10  MHz  is 
that  it  is  intractable  to  use  a  bandlimiting  FIR  filter  response  that  accurately  restricts 
the  message  bandwidth  to  +/{,  Hz.  A  consequence  of  this  is  that  although  representative 
relative  accuracy  between  OOK  and  FSK  is  shown  in  Figure  4  there  will  be  some  absolute 
error.  Nevertheless,  comparison  of  the  OOK  TDOA  accuracy  at  8.5  dB  in  Figure  4  and  0 
dB  in  Figure  3  reveals  fairly  close  consistency. 

Figure  4  is  depicted  with  a  log  vertical  scale  to  better  accommodate  the  full  range  of 
TDOA  estimation  accuracies  encountered.  It  shows  a  fairly  constant  ratio  between  the 
mis  TDOA  accuracies  for  OOK  and  FSK  estimation  at  high  SNR,  with  the  accuracy  im¬ 
proving  with  increasing  frequency  deviation  /o-  However  at  low  SNR,  an  almost  reverse 
relationship  is  evident.  Both  of  these  observations  are  intuitively  satisfying;  at  high  SNR 
one  expects  the  accuracy  to  improve  with  increasing  frequency  deviation  due  to  the  sharp¬ 
ening  of  the  correlation  function,  whereas  the  associated  ambiguities  in  the  correlation 
function  are  closer  in  both  time  and  magnitude  to  the  correct  peak,  leading  to  poorer  low 
SNR  performance. 

Prior  to  considering  the  quantitative  comparison  between  OOK  and  FSK  performance, 
it  should  be  first  considered  that  the  rms  TDOA  shown  in  Figure  4  is  plotted  against  SNR 
within  the  signal  bandwidth,  which  is  approximately  10  kHz  for  the  OOK  signal  and 
twice  this  for  FSK.  Computed  at  +8.5  dB  SNR,  the  reduction  in  rms  TDOA  achieved 
in  simulation  through  using  FSK  in  lieu  of  OOK  is  8.4,  15.2,  27.3  and  80.0,  for  FSK 

6That  is,  spurious  peaks  located  some  distance  from  the  correct  peak  are  sometimes  encountered  in  one 
or  more  simulation  runs  out  of  100. 

'  Equation  (4)  shows  that  for  the  case  where  the  signal  and  noise  are  white,  a  reduction  in  L  can  be  offset 
by  a  corresponding  increase  in  equivalent  SNR  7.  Application  of  equation  (5)  allows  the  corresponding 
SNR  to  be  evaluated.  Specifically,  if  the  SNR  for  each  received  signal  is  7  =  0  dB  then  7 eq  =  —1.76  dB. 
A  10  dB  increase  in  equivalent  SNR  is  7 eq  =  8.24  dB.  which  corresponds  to  7  «  8.5  dB.  It  is  conjectured 
that  although  this  result  is  only  specifically  applicable  to  white  data,  it  should  also  be  suggestive  for  the 
signal  of  interest. 
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Figure  1:  Magnitude  complex  correlation  function,  OOK  and  FSK.  fs  =  312.5kHz 
SNR=mdB,  fb  =  5  kHz. 


x  104 


Fiqure  2:  Maqnitude  complex  correlation  function,  OOK  and  FSK.  =  312.5kHz 
SNR=—9dB,  fh  =  5  kHz. 
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frequency  deviation  equal  to  5,  10,  20  and  40  kHz  respectively.  For  the  current  simulations 
with  fb  =  5  kHz,  B  =  10  kHz,  evaluation  of  equation  (27)  suggests  a  relative  accuracy 
equal  to  approximately  ctt,ook / &t,2FSK  =  2, 3.6,  7  and  13.9  respectively.  Quantitative 
consistency  between  the  simulation  results  and  theory  is  not  clearly  evident,  although  the 
trends  are  common.  The  reason  for  this  inconsistency  is  due  to  the  assumption  made 
in  the  theoretical  development,  that  the  signal  has  a  bandlimited  white  power  spectrum. 
Just  as  was  identified  and  resolved  in  section  2.2  for  the  study  of  real,  passband  TDOA 
estimation,  the  degree  of  improvement  in  TDOA  accuracy  is  improved  substantially  when 
the  signal  has  a  non-white  power  spectrum. 

Both  the  approximate  theory  and  simulation  results  show  a  potentially  substantial 
improvement  to  the  rms  TDOA  estimation  accuracy  through  using  an  FSK  modulation 
in  preference  to  OOK.  For  the  simulation  case  of  Figure  4,  when  the  SNR  exceeds  around 
+5dB,  a  significant  improvement  is  evident  for  all  frequency  deviations  considered.  How¬ 
ever,  if  it  is  desirable  to  ensure  a  performance  improvement  across  all  SNR  scenarios  exam¬ 
ined  in  the  Figure  4  simulations  (that  is,  down  to  as  low  as  -3.5  dB),  an  FSK  modulation 
with  /o  =  fb  =  5  kHz  frequency  deviation  can  be  used. 

The  results  from  Figure  4  can  be  used  to  infer  the  performance  of  the  original  scenario 
of  interest,  with  an  integration  period  of  1.5  seconds  by  translating  the  horizontal  axis  by 
8.5  dB.  That  is,  the  performance  for  the  1.5  second  scenario  at  0  dB  SNR  should  be  equal 
to  that  in  Figure  4  at  +8.5  dB.  For  greater  SNR  levels,  performance  equivalence  will  be 
achieved  at  an  SNR  difference  tending  to  10  dB,  but  at  lower  SNR,  the  difference  will  be 
less8. 


4  Conclusion 

This  report  has  explored  the  potential  for  achieving  relative  TDOA  accuracy  improvement 
through  using  FSK  signalling  in  lieu  of  OOK.  It  has  been  shown  through  both  analysis  and 
simulation  that  a  significant  improvement  in  accuracy  can  be  achieved  at  an  appropriately 
high  SNR.  The  improvement  in  the  rms  TDOA  error  can  exceed  an  order  of  magnitude 
for  signals  and  conditions  of  practical  interest. 
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Figure  3:  rms  TDOA  for  OOK  and  FSK  modulation.  fs  =  312.5 kHz  and  f(  =  5  kHz. 
Note:  SNR  is  measured  within  signal  bandwidth  (for  FSK  this  is  ±/&  bandwidths  around 
each  of  fc  ±  f0) 


SNR  (dB)  in  approximate  signal  bandwidth 

Figure  f:  rms  TDOA  for  OOK  and  FSK  modulation.  fs  =  10 MHz,  f),  =  5  kHz  and 
L  =  0.15  sec.  Note:  SNR  is  measured  within  approximate  signal  bandwidth  (for  FSK ,  this 
is  ±fb  bandwidths  around  each  of  fc  ±  ,/q ) 
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Appendix  A  Alternative  exploitation  of 

correlation  phase 

A.l  Improved  Complex  Baseband  TDOA  Estimation 

The  TDOA  estimation  benefits  of  real  passband  processing  can  be  also  captured  in  the 
complex  baseband  domain  by  exploiting  the  phase  information  within  the  complex  corre¬ 
lation  function.  There  are  two  approaches  to  this: 

Approach  1  Define 

R{t)  =  R(t)  exp(— j27t(/c  -  /j)t)  (Al) 

=  exp(j'27r/c(T  -  r))  J s(t  -  r)s*(t  -  T)dt,  (A2) 

where  R(t)  is  defined  in  equation  (11). 

Estimation  of  the  TDOA  is  then  via: 

T  =  argrnax  9?[.R(t)] 

T 

=  argrnax  [cos(27 r/c(T  —  r))  J ' a(t  —  r)a(t  —  T)  +  b(t  —  r)b(t  —  T)dt 
+  sin(27r/c(T  —  r))  ~  T)b{t  —  T)  —  b(t  —  r)a(t  —  T)dt ]. 


Approach  2  This  approach  involves  two  steps: 

(i)  estimate  TDOA  from  the  magnitude  of  the  complex-baseband  correlation  function: 


T\  =  argrnax  \R(t)  . 

T 

(A3) 

(ii)  Define: 

=  phase  [i?(Ti )] 

(A4) 

Te  =  T  —  T\. 

(A5) 

Then 

(j)\  =  27 r(/cT  —  f [T i )  +  n2n,  where  n  =  0,-1, ... 

(A6) 

=  27 r((/c  -  fj )7"l  +  fcTe )  +  77.27T. 

(A7) 

Hence,  estimation  of  Te  may  be  obtained  via: 

fe  =  (</>1/(2n)-n)-(fc-fI)f1)/fc 

(AS) 

leading  to  the  final  TDOA  estimate: 

f  =  Ti  +  fe. 

(A9) 
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It  should  be  acknowledged  that  both  of  these  approaches  are  susceptible  to  the  n//c 
estimation  ambiguities  also  inherent  to  the  real  passband  processing  approach.  It  should 
also  be  pointed  out  that  both  of  these  complex  baseband  approaches  above,  and  the  real 
IF  approach  in  Appendix  A. 2,  require  a  sampling  rate  fs  >  2  fc  to  obtain  the  phase  related 
TDOA  estimation  improvements. 


A. 2  Real  IF  TDOA  Estimation 


The  improved  TDOA  estimation  accuracy  obtainable  from  the  (real  and  complex  process¬ 
ing)  methods  in  Appendix  A.l  can  be  also  obtained  with  processing  the  real  signals  down 
converted  to  a  suitable  intermediate  frequency  (IF). 

u\(t) — >  x\(t)  =  a(t)  cos(2tt fjt)  —  b(t)  sm(2n fit) 

U2(t) — >  x2(t)  =  a(t  —  T)  cos(27r(//t  —  fcT))  —  b(t  —  T)  sin(27r(//f  —  fcT)). 


Define: 


xi  (t  -  t) 

where  x\ n{t) 
R(r) 


a(t  -  t)  cos(27r(//f  -  /cr))  -  b(t  -  r)  sin(27r(//t  -  fcr)) 
x\ (t  -  r)  cos(27 r(/c  -  //)r)  +  x1H{t  -  r)  sin(27r(/c  -  //)r) 
Hilbert  transform[x'i  (t)] 

J  x±(t  —  r)x2(t)dt. 


Improved  TDOA  estimation  is  then  obtained  via: 

T  =  argrnax  R(t) 

T 

=  argrnax  [cos(2tt fc(T  —  t))  J a(t  —  r)a(t  —  T)  +  b(t  —  r)b(t  —  T)dt 
+  sin(27r/c(T  —  r))  ~  T)b(t  —  T)  —  b(t  —  r)a{t  —  T)dt]. 
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